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Abstract. We have carried out a molecular survey of the Class IM protostar NGC 7129 - FIRS 2 (hereafter 
FIRS 2) and the Herbig Be star LkHa 234 with the aim of studying the chemical evolution of the envelopes of 
intermediate-mass (IM) young stellar objects (YSOs). Both objects have similar luminosities (~500 Lq) and are 
located in the same molecular cloud which minimizes the chemical differences due to different stellar masses or 
initial cloud conditions. Moreover, since they are located at the same distance, we have the same spatial resolution 
in both objects. A total of 17 molecular species (including rarer isotopes) have been observed in both objects and 
the structure of their envelopes and outflows is determined with unprecedent detail. 

Our results show that the protostellar envelopes are dispersed and warmed up during the evolution to become a 
pre-main sequence star. In fact, the envelope mass decreases by a factor >5 from FIRS 2 to LkHa 234, while the 
kinetic temperature increases from ~13 K to 28 K. On the other hand, there is no molecular outflow associated 
with LkHa 234. Fhe molecular outflow seems to stop before the star becomes visible. 

These physical changes strongly affect the chemistry of their envelopes. The N2H + and NH3 abundances seem to 
be quite similar in both objects. However, the H 13 CO + abundance is a factor of ~3 lower in the densest part of 
FIRS 2 than in LkHa 234, very likely because of depletion. In contrast, the SiO abundance is larger by a factor of 
~100 in FIRS 2 than in LkHa 234. CS presents a complex behavior since its emission arises in different envelope 
components (outflow, cold envelope, hot core) and could also suffer from depletion. The CH3OH and H2CO 
column densities are very similar in FIRS 2 and LkHa 234 which implies that the beam averaged abundances are 
a factor >5 larger in LkHa 234 than in FIRS 2. The same case is found for the PDR tracers CN and HCN which 
have similar column densities in both objects. Finally, a complex behavior is found for the deuterated compounds. 
While the DCO + /H 13 CO+ ratio decreases by a factor of ~4 from FIRS 2 to LkHa 234, the D 2 CO/H 2 CO ratios is 
within a factor 1.5 in both objects. The detection of a warn CH3CN component with Tfe>63 K shows the existence 
of a hot core in FIRS 2. Thus far, only a handful of hot cores have been detected in low and intermediate mass 
stars. 

Based on our results in FIRS 2 and LkHa 234, we propose some abundance ratios that can be used as chemical 
clocks for the envelopes of IM YSOs. The SiO/CS, CN/N 2 H+, HCN/N 2 H+, DCO+/HCO+ and D 2 CO/DCO+ 
ratios are good diagnostics of the protostellar evolutionary stage. 

Key words. Stars: formation - pre-main sequence - individual: LkHa 234 - ISM: abundances - clouds - individual: 
NGC 7129 



1. Introduction 

Chemistry is a powerful tool to study young stellar objects 
(YSOs) and their environments. On the one hand, chem- 
istry is a diagnostic tool of the different envelope com- 
ponents. On the other hand, it is a good time indicator 
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during the protostellar evolution. Chemical studies have 
been used t o determine the physical structure of low- mas s 
YSOs (e.g. lMaret et al. 1120041 [lortrensen et al. Il2004albl h 

These studies have revealed, for example, the presence of 
warm regions where the ices evaporate giving rise to re- 
gions similar to hot cores in massive protostars. However, 
the chemistry of the two classes of objects is different 
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Fig. 1. Map of the 1.3mm continuum flux toward NGC 7129. Levels are 22.3, 44.6, 89.2, 178.4 to 713.6 mJy/beam by 
178.4 mJy/beam. Crosses indicate the millimeter sources and stars the infrared sources. The red and blue contours 
represent the redshifted and blucshiftcd high-velocity gas of the nebula as traced by the 12 CO J=2^1 line. Velocity 
intervals are [-30,-13] kms -1 for the blue lobe (grey contours) and [-7,11] kms -1 for the red-lobe (dark contours). 
Contours are 10 to 60 K kms -1 in steps of 10 K kms -1 for the blue lobe and 10 to 100 K kms -1 in steps of 15 K 
kms -1 for the red lobe. 



and rises questions on the mechanisms that lead to the 
observed che mical complexity and its dependence on th e 
stellar mass (|Cazaux et al. Il2003t Eottinelli et al. 112004ft ■ 

Chemistry has also been used as a time evolution 
indicator in both low-mass a nd high-mass objects. For 
instance, iMaret et al. I l|2004f) studied a sample of low- 
mass Class objects and found some indications that the 
H2CO abundance increases during the protostellar evolu- 
tion. However, the different abundances can be due to dif- 
ferent initial cloud conditions from which the protostars 
evolved. Further chemical studies are required to clarify 
this subject. 

In this paper, we present a chemical study of the two 
intermediate mass (IM) YSOs, FIRS 2 and LkHa 234. 
Contrary to low-mass protostars, IM YSOs have been very 
little studied thus far, specially young proto stars. For ex- 
ample , out of 42 Class sources compiled bv lAndre et al. I 
(2000), only six had luminosities in excess of 40 L Q (the 
precursors to HAe stars) and only one had a luminosity 
of ~ 10 3 L (the precursor to HBe stars). IM YSOs (M* 
~ 2 -10 M©) are not only an important link between low 



mass and high mass stars but they share many character- 
isitcs of high mass star formation (clustering, PDR/HII 
regions) without the disadvantage of being too distant (d 
~ 1 Kpc), or too complex. They are also important for 
the understanding of planet formation since Herbig Ae 
stars are the precursors of Vega-type systems. On a larger 
scale, th ey dominate the mea n UV interstellar field in our 
Galaxy llWolfire et al. 112003ft . 

2. Target selection and observational strategy 

FIR S 2 has been classified as a Class IM object 
l|Eiroa et al. 1 1199^1 and is, very likely, the youngest IM 
object known at present. An energetic bipolar m olecular 
outflow is associated with it IjFuente et al. I 2001). LkHa 
234 is a embedded HBe star which still keeps a massive 
envelope (M~16 M Q ) b ut no bipolar molecular outflow 
seems to be driven by it l)Fuente et al. ll200lU2002ft . These 
objects have the peculiarity of having similar luminosities 
(~ 500 L Q ) and be located in the same molecular cloud. 
During the protostellar and pre-main sequence evolution, 
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the lumin osity remains quite constant for a given stel- 
lar mass llAndre et al. 2000). Thus, the same luminos- 
ity implies similar stellar mass. In addition, both sources 
are located in the same molecular cloud. This minimizes 
chemical effects due to very different stellar mass (both of 
them have M* ~ 5 M Q ) and/or different initial cloud con- 
ditions. Thus, in spite of our reduced sample, we have an 
excellent opportunity to find an evolutionary track for IM 
YSOs. In Fig. 1 we show the 1.3mm continuum map of the 
reflection nebula NGC 7129. Both, FIRS 2 (NGC 7129- 
FIRS 2) and LkHa 234 (NGC 7129 - FIRS 1) are associ- 
ated with intense centrally peaked 1.3mm continuum and 
far infrared sources. We have also overlaid the red and blue 
lobes (as traced by the 12 CO J=2^1 line) of the bipolar 
outflows associated to these YSOs (see Fuente et al. 2001 
for a more detailed description). Clust ering becomes sig - 
nificant in this range of stellar masses ijTesti et al. l fl999) . 
Interferometric observations towards LkHcv 234 show the 
existence of a young infrared companion (IRS 6) which 
is very likely the exciting source of th e bipolar m o lecula r 
outflow and the [SII] jet detected bv iRav et al. I l|l990l) . 
The quadrupolar morphology of the ouflow detected in 
NGC 7129-FIRS 2 is also due to the superposition of two 
bipolar molecular ou t flows FIRS 2-out 1 an d FIRS 2-out 2 
I Fuente et ai~ll200ll) . iFuente et all (|200lh proposed that 
FIRS 2-out 1 is associated with the Class protostar while 
FIRS 2-out 2 is more likely associated with a more evolved 
infrared star (FIRS 2 - IR). 

A complex chemical evolution occured in the YSO en- 
velopes during the protostellar evolution. This involves 
accretion of species in an icy mantle during the pre- 
collapse phase, followed by grain-surface chemistry and 
evaporation of ices on ce the YSO has st arted to heat 
its surroundings (e.g. iBrown et al. I ^988). In massive 
YSOs, the evaporated molecules drive a rapid high- 
temperature chemistry for a period of 10 4 -10 5 years, 
resulting in the complex, saturated organic molecules 
(CH 3 OCH 3 ,CH3CN, C?HbOH,...), which are characteris- 
tic of a hot core (e.g. ICharnlev et al. Ill902i ICaselli et al. I 
19931 iRodeers fc Charnlevll2003t iNomura fc Millar Il2004i 
Viti et alTI l2004h . Once most of the envelope has been 
cleared up, the UV radiation can escape to form a photon- 
dominated-region (PDR) and, in the case of massive 
stars, an HII region. Simultaneously, the energetic bipolar 
outflows develope a shock chemistry in the surrounding 
molecular cloud. 

To discern between the different envelope components 
and determine the protostellar envelope evolution, we 
have selected a set of molecular tracers. Specifically, NH3, 
N 2 H+, H 13 CO+ and HC 18 0+ have been observed as 
tracers of the cold envelope; the volatile species CH3OH 
and H2CO to study the warm envelope where the icy 
grain mantles have been evaporated; the complex molecule 
CH3CN to trace the hot core; and, CN and HCN in or- 
der to trace the incipient PDR. In addition, we have also 
observed several CS and C 34 S lines that are useful to con- 
strain the physical conditions of the envelope and SiO as 
an excellent tracer of the shock chemistry associated to 



the outflow. The observation of the deuterated compounds 
N2D + ,DCO + and D2CO also provide us some informa- 
tion about how the deuterium fractionation is affected 
by the YSO evolution. Obviously, the "molecular tracer" - 
"envelope component" correspondence is not unique and 
all the species have contributions from other envelope 
components. Arguments based on the morphology of the 
emission, kinematics and different excitation conditions 
are used to discern between the different components in 
these cases. 

3. Observations 

The (J,K) = (1,1), (2,2), (3,3) and (4,4) inversion lines 
of ammonia were observed using the Effelsberg 100-m ra- 
diotelescope of the MPIfR in December 2000 and October 
2002. We observed only one position in both sources. The 
half power beam width (HPBW) of the telescope at the 
rest frequency of the NH 3 lines, 23.7 GHz, was 40". We 
used the new cooled dual-channel HEMT K-band receiver 
with a typical system temperature of 200 K on a main 
beam brightness temperature scale. The 8192-channel au- 
tocorrelator was used as the backend. The four ammonia 
lines were observed simultaneously with a total bandwidth 
of 10 MHz and a channel separation of 0.098 kms -1 . We 
estimate that line intensities are accurate to within ±15%. 

The main set of observations were carried out with 
the IRAM 30m telescope in Pico de Veleta (Spain) during 
three different observing periods in May 1999, July 2002 
and August 2003. The list of observed lines, the telescope 
characteristics at each frequency and a summary of the 
observations are shown in Table 1. When possible, all the 
lines of the same molecule were observed simultaneously 
in order to avoid observational errors. The backends were 
an autocorrelator split in several parts and a 256 x 100 
kHz filter-bank. All the lines were observed with a spec- 
tral resolution of ~ 78 kHz except in the cases that are 
explicitly indicated in Table 1. The intensity scale used in 
this paper is main brightness temperature. Comparing the 
intensity of some pattern lines in different observing peri- 
ods, we estimate that the calibration is accurate within a 
20 % at 3mm and within a 30% at 1.3mm. We have ob- 
served a high S/N ratio spectra towards the star positions 
for all the lines listed in Table 1 . In addition, we have car- 
ried out small maps in the most intense lines (see Table 
1). In the case of the SiO 2->l, SiO 3^2, H 13 CO+ 1^0 
and CS 3— >2 lines we have mapped the entire outflows 
as traced by the high velocity CO emission. Fits to the 
oberved lines at the (0,0) position are shown in Tables 2 
and 3. 

4. Analysis 

The data have been analysed using the rotation diagram 
method. The molecular constants, the upper state en- 
ergies and the partition functions required for applying 
this method have been taken from the JPL line catalog 
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Table 1. Description of the IRAM 30m observations 



Line 


Freq. (GHz) 


HPBW 


T]MB 


Observed positions 


N 2 H + l->0 


93173.2 


26.5" 


0.71 


60"x60" map 


CN l-»0 


113490.0 


23" 


0.65 


60"x60" map (FIRS 2), 60" x 100" map (LkHa 234) 


HCN l->0 


88631.8 


27.5" 


0.73 


60"x60" map (FIRS 2), 60" x 100" map (LkHa 234) 


H 13 CO+ l->0 


86754.3 


27.5" 


0.73 


complete map 


HC 1S + l->0 


85162.2 


27.5" 


0.73 


(0,0) 


SiO 2-s-l 


86846.9 


27.5" 


0.73 


complete map 


SiO 3^2 


130268.6 


20" 


0.60 


complete map 


CS 3^2 


146969.0 


16.5" 


0.55 


complete map 


C 34 S 2->l 


96412.9 


25.9" 


0.70 


(0,0) 


C 34 S 3^2 


144617.1 


17" 


0.55 


(0,0) 


C 34 S 5^4 


192818.5 


13" 


0.47 


(0,0) 


CH 3 OH 2(1,3)^1(1,3) 


96755.5 


26" 


0.70 


72"x72" map 


CH3OH 2(0,3)^1(0,3) 


96744.6 


26" 


0.70 


72"x72" map 


CH3OH 2(1,4)^1(1,4) 


96741.4 


26" 


0.70 


72"x72" map 


CH3OH 2(0.1)^1(0.1) 


96739.4 


26" 


0.70 


72"x72" map 


CH3OH 3(0,3)^2(0,3) 


145093.7 


17" 


0.55 


(0,0) 


CH3OH 3(1,4)->2(1,4) 


145097.5 


17" 


0.55 


(0,0) 


CH3OH 3(0,1)^2(0,1) 


145103.2 


17" 


0.55 


(0,0) 


CH3OH 3(2, 2)->2(2, 2) a 


145124.4 


17" 


0.55 


(0,0) 


CH3OH 3(1,3)^2(1,3) 


145131.9 


17" 


0.55 


(0,0) 


CH3OH 5(2,4)^4(2,4) 


241904.1 


10" 


0.40 


72"x72" map 


CH3OH 5(0.1)^4(0.1) 


241791.4 


10" 


0.40 


72"x72" map 


CH3OH 5(1,4)^4(1,4) 


241767.2 


10" 


0.40 


72"x72" map 


CH3OH 5(3,1)^4(3,1) 


241832.9" 


10" 


0.40 


72"x72" map 


CH3OH 5(2,2)^4(2,2) 


241842.3 


10" 


0.40 


72"x72" map 


CH3OH 8(1,4)^7(0,3) 


229758.7 


10.5" 


0.42 


(0,0) 


CH3OH 11(1,4)^10(2,4) 


104300.5 


24.5" 


0.68 


(0,0) 


H2CO 3i2 — ^2n 


225697.8 


10.5" 


0.43 


60"x60" map (FIRS 2), 72"x72" map (LkHa 234) 


H2 CO 2 12 — ^ In 


140839.5 


18" 


0.57 


60"x60" map (FIRS 2), 72"x72" map (LkHa 234) 


H2^^CO 2i2 — ^ln 


137449.9 


18" 


0.57 


(0,0) 


CH 3 CN 5(0)-»4(0) 


91987.0 


27" 


0.72 


(0,0) 


CH 3 CN 5(1)^4(1) 


91985.3 


27" 


0.72 


(0,0) 


CH3CN 5(2)^4(2) 


91980.0 


27" 


0.72 


(0,0) 


CH 3 CN 5(3)^4(3) 


91971.4 


27" 


0.72 


(0,0) 


CN 2-t-l 


226874.7 


11" 


0.42 


60"x60" map (FIRS 2), 60" x 100" (LkHa 234) 


N 2 D+ 3^2 


231321.7 


10.5" 


0.42 


(0,0) 


DCO+ 2->l 


144077.3 


17" 


0.56 


60"x60" map 


DCO+ 3^2 


216112.6 


11" 


0.44 


60"x60" map 


D 2 CO 4o4^3 03 


231410.3 


10.5" 


0.42 


(0,0) 



a Only observed with a frequency resolution of 1 MHz. 



(Pic kett et al. 1998). This method gives the rotation tem- 
perature and total column density of a particular species 
knowing the integrated line intensities of several lines with 
different upper state energies. The rotation temperatures 
and column densitites estimated in this way are shown in 
Table 4. The rotation temperature is a lower limit to the 
gas kinetic temperature, and only constitutes a good mea- 
sure of the kinetic temperature if the lines are thermalized. 
This is the case of the ammonia inversion lines which are 
thermalized with densities, n^lO 3 cm~ 3 . For this reason, 
we have used the NH3 inversion lines to estimate the gas 
kinetic temperature of the different envelope components. 

For some molecules we have only observed one tran- 
sition. In this case we have calculated the total column 
density assuming optically thin emission, local thermody- 
namic equilibrium (LTE) and the rotation temperature 



derived from a molecule with similar excitation require- 
ments. In these cases we have marked with the superindex 
"a" the rotation temperature in Table 4. 

The molecules NH 3 , N 2 H+, CN and HCN present hy- 
perfine splitting. This allows us to derive the line opacity 
directly from the hypcrfmc line ratios. In these cases the 
column densities have been estimated directly from the 
line opacities. For N2H + and NH3, we have also calcu- 
lated the source size assuming the rotation temperature 
derived from the NH 3 inversion lines when the opacity of 
the main component is determined. We have made LVG 
calculation for SiO, CS and C 34 S. These calculations are 
shown in Tables 6 and 8. In all the cases, we have fitted 
the densities assuming a fixed kinetic temperature. The 
assumed gas kinetic temperatures are based on those de- 
rived from the NH3 data in Section 5.1 and are shown in 
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Line 




V (km s" 1 ) 


AV (km s- 1 ) 




NH 3 (1,1) 


3.08(0.07) 


-9.59(0.01) 


1.0(0.1) 


0.1 a 


NH 3 (2,2) 


0.70(0.02) 


-9.63(0.02) 


1.3(0.1) 


0.1 a 


NH 3 (3,3) 


0.12(0.01) 


-9.35(0.18) 


3.8(0.5) 


0.1 a 


NH 3 (4,4) 


rms — 


0.05 K kms" 1 


with Av— 3.0 kms" 


-1 


N 2 H+ 1-0 


9.22(0.02) 


-9.51(0.01) 


1.91(0.01) 


1.28(0.10) 


CN 1-0 


2.3(0.3) 


-9.7(0.1) 


1.4(0.1) 


1.9 


HCN 1—0 


2.3(0.3) 


-9.3(0.1) 


3.1(0.1) 


0.7 


Line 


Area (K kms™ 1 ) 


V (km s" 1 ) 


AV (km s" 1 ) 


T M5 (K) 


H la CO + 1—0 


1.92(0.02) 


-9.74(0.01) 


1.27(0.02) 


1.39 


HC ls O + 1—0 


0.17(0.02) 


-9.68(0.09) 


1.1(0.2) 


0.14 


SiO 2 — 1 


1.55(0.09) 


-7.5(0.2) 


6.8(0.4) 


0.21 




0.23(0.06) 


-9.5(0.1) 


1.6(0.4) 


0.13 


SiO 3—2 


2.50(0.20) 


-7.3(0.2) 


6.9(0.7) 


0.34 


CS 3—2 


3 89(0 1 2^ 


-9.1(0.1) 


5 81 (0 ?4'i 


0.60 




1.88(0.09) 


-9.9(0.1) 


1.19(0.04) 


1.48 


n 1 fifn 

U.1U l^U.UO ) 


-7 8(0 8) 


3 1(1 3) 


05 


O o o — >z 


n 9S/TI 07 s ! 


-9 8(0 1) 


1 3(0 2) 


20 


l.Uo(U.Uo ) 


-9.4(0. 1) 


4.3(0.4) 


0.23 


C 34 S 5 >4 


1.21(0. lb) 


-8.8(0.3) 


4.2(0.7) 


0.27 


CH3UH 2(1, 3) — 3) 


0.6(0.3) 


-9(1) 


5(3) 


0.11 


Orl3Url 2(0,3)-— >1(0, 3) 


1.85(0. 0b ) 


-8.20(0.09) 


6. 1(0.2) 


0.28 


CH3UH 2(1,4) — ► 1(1,4) 


4.71(0.05) 


-9.2(0.2) 


3.0(0.2) 


1.47 


CH3OH 2(0,1)— 9-1(0, 1J 


5.83(0.07) 


-9.58(0.01) 


5.3(0.1) 


1.03 


CH3OH 3(0,3)— »2(0, 3) 


4.12 (0.04) 


-8.78(0.08) 


6.98(0.08) 


0.55 


CH3OH 3(1,4)^2(1,4) 


8.9(0.09) 


-8.32(0.01) 


5.71(0.07) 


1.46 


Orl3Url 3(0, 1)— »2(0, 1) 


9.6(0. 1) 


-8.34(0.02 ) 


5.07(0.06) 


1.78 


CH3OH 3(2,2)— >2(2, 2) 


2(1) 


-12(4) 


7(10) 


0.27 


Orl3Url 3(1,3) — >2(1,3) 


1.9(0. 1) 


-8.7(0.2 ) 


6.4(0.5) 


0.28 


L/H3UM 5(2,4) — >4(2,4) 


3(1) 


-8.9(0.7; 


a /n't 

4(2) 


0.70 


riTT nu c/n 1 '\ a frt 1 \ 
L/H3UM 5(0, 1)— >4(0, 1) 


12.3(0.7; 


-7.4(0. 1) 


4.6(0.3) 


2.51 


PTT OFF rr/i y] \ ./l/l /l^ 

O Fl 3 FJ H 0(1,4: ) — »4( 1 ,4: J 


O.0(U.O ) 


•7 { r\ n\ 
- 1 .o(U.z ) 


a afn k\ 
4.D(U.OJ 


1 .80 


piJ OFT tifQ 1^ ./l/ 1 } 1^ 






O-'J(I) 


U . z 


pij nw ^fo o^-^/i/^) 9^ 

V_ Fl 3 rl uyZ^Z ) - — >4 l^z , z j 


±.o( 1 .0) 


-10(4_) 


((llj 


n 17 
u. 1 / 


OFF OFF fifl /I .Tl'n 

FI3 FJFL 0(1,4:) — '/(U,Oj 


4. (U.4J 


-7.7(0.1) 


4.3(0.3) 


0.98 


OFF OFF 1 1 h /I \ jn/O A\ 

OII3 Wll 11(1,4 ) — >■ 10(^,4, 


U. 14(11.114 J 


! n a \ 
-y.y(u.4 ) 


3.6(5.9) 


0.04 


H2L'U 012^^-^11 


9 ^fn 9*1 
z .0 . z ) 


-9 6(0 1) 


1 3(0 1) 


1 66 




6.7(0.3) 


-8.6(0.2) 


7.7(0.5) 


0.81 


H9CO 2l9— >ll 1 


3.0(0.1) 


-9.66(0.01) 


1.41(0.03) 


2.00 




6.0(0.1) 


-8.97(0.06) 


6.3(0.2) 


0.91 


H2 13 CO 2i2 — >ln 


0.33(0.04) 


-9.4(0.2) 


2.8(0.5) 


0.11 


CH 3 CN 5(0)— 4(0) 


0.24(0.05) 


-9.1(0.4) 


4.2(0.8) 


0.05 


CH 3 CN 5(1)— 4(1) 


0.25(0.05) 


-9.1(0.6) 


5(1) 


0.05 


CH 3 CN 5(2)— 4(2) 


0.11(0.02) 


-9.7(0.9) 


4(2) 


0.02 


CH 3 CN 5(3)— 4(3) 


0.12(0.02) 


-8(1) 


7(3) 


0.02 


CN 2—1 


2.74(0.24) 


-9.7(0.1) 


2.6(0.2) 


1.0 


N 2 D+ 3—2 


0.64(0.05) 


-9.90(0.04) 


1.1(0.1) 


0.53 


DCO+ 2 — 1 


2.93(0.04) 


-9.86(0.01) 


1.0(0.1) 


2.61 


DCO + 3—2 


1.9(0.1) 


-9.73(0.02) 


1.0(0.1) 


1.83 


D 2 CO 4 04 - 3 03 


0.3(0.1) 


-9.6(0.3) 


3.9(1) 


0.05 



a In case of optically thin emission, T m cannot be determined and is set arbitrarily to 0.1. 



the (4,4) line. We have fitted these lines using the NH3 
procedure in the CLASS (this is the program dedicated 
to the analysis of single-dish spectra of the GILDAS 
softwares ( |http: / / www.iram.fr /IRAMFR/GILDAS ) . This 
procedure fit all the hyperfine components assuming equal 
excitation temperature, central velocity and linewidth. 
The parameters given by the procedure are (T ex — 
Tb g )xT m , V, AV, and r m where r m is the main group 
opacity and T ex the excitation temperature (see Bachillcr 
et al. 1987 for a more detailed description of this pro- 
cedure). In case of optically thin emission r m cannot be 
determined and is set arbitrarily to 0.1. 

The fits to the NH 3 lines in NGC 7129-FIRS 2 are 
shown in Table 2. The NH3 emission is optically thin in 
all the lines. However, the linewidth of the (3,3) line is 
almost a factor of 3 larger than the linewidth of the (1,1) 
and (2,2) lines (see Table 2 and Fig. 2). In addition, the 
central velocity of the lines is slightly lower than those of 
the lower energy lines. This suggests that the (3,3) line is 



Table 6 and 8. We have used the CS col lisional coefficients 
calculated bv lGreen k. Chapmanl ljl978|) in the LVG calcu- 
lations. The same collisional coefficients are used for SiO. 
This is a reasonable approximation since both molecules 
have the same mass and similar dipole moments. In the 
case of CS, we have been able to calculate the opacity and 
the source size because we have observed the main iso- 
tope and the rarer isotope C 34 S. All the column densities 
in Table 4, 6 and 8 are beam averaged column densities. 

5. Results 

5.1. NH 3 

We have observed the ammonia (1,1), (2,2), (3,3) and 
(4,4) inversion transitions towards NGC 7129 - FIRS 2 
and LkHa 234. All the spectra are shown in Fig. 1. The 
(1,1), (2,2) and (3,3) lines have been detected in both 
sources but only an upper limit has been obtained for 
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Table 3. Observational parameters towards LkHa 234 



Line 


TMB XT m 


V (km s" 1 ) 


AV (km s- 1 ) 




NH 3 (1,1) 


0.60(0.03) 


-9.87(0.02) 


0.9(0.2) 


0.86(0.07) 


NH 3 (2,2) 


0.27(0.04) 


-9.94(0.05) 


1.1(0.2) 


0.1 Q 


NH 3 (3,3) 


0.10(0.03) 


-11(1) 


1.6(0.7) 


0.1 a 


NH 3 (4,4) 


rms 


= 0.13 K kms -1 with Av=l.l kms" 


1 


N 2 H+ 1—0 


1.53(0.06) 


-9.48(0.01) 


1.82(0.04) 


1.32(0.7) 


CN 1—0 


1.9(0.2) 


-9.9(0.1) 


2.6(0.1) 


0.8 


HCN 1—0 


3.2(0.3) 


-10.0(0.1) 


2.9(0.1) 


0.1 a 


Line 


Area (K kms" 


" 1 ) V (km s -1 ) 


AV (km s _I ) 


T MB (K) 


H la CO + 1—0 


1.00(0.03) 


-y.oui u.uo j 


1 . 1 y ( u . u 1 ) 


0.52 


HC ls O + 1—0 


0.14(0.04) 


-10. 1(0. 2 j 


2(1) 


0.06 


SiO 2 — 1 


0.31(0.06) 


-9.9(0.9) 


8.5(2.0) 


0.03 


giO 3—2 


0.60(0.08) 


-8.7(0.2 ) 


3.2(0.5) 


0.18 


CS 3—2 


14.9(0.2) 


-9.7(0.1) 


3.0(0.1) 


4.67 


C 34 S 2—1 


66f0 04 s ! 


-9.9(0.1) 


3.4(0.2) 


18 


p34q q .n 


1 A 1 t C\ C\K\ 
1.4HU.UO ) 


-9.5(0.1) 


3.2(0.2) 


0.41 


C 34 S 5—4 


1.71(0.22) 


-9.1(0.3) 


4.5(0.8) 


0.36 


CH3CJII 2(1, 3J — *l{l,o) 


0.3(0.1) 


-7(3) 


12(5) 


0.02 


OH3UH 2(u,3) — »l(u,3) 


0.2u(0.04) 


-9.8(0.3) 


2.9(0.7) 


0.06 


CH3UII 2(1,4)— >1(1, 4) 


0.48(0.03) 


-9.6(0.1) 


2.2(0.2) 


0.20 


Orl3Url z(U, 1J — >1(U,1) 


0.64(0.04) 


-9.5(0.1) 


5.2(0.4) 


0. 11 


CH3OH 3(0,3)— »2(0, 3) 


1.27(0.09) 


-9.5(0.1) 


4.6(0.4) 


0.26 


CH3CJII 3(1,4)— >2(1, 4) 


1.91(0.09) 


-9.2(0.1) 


4.8(0.3) 


0.37 


OH3UH 3(u, 1)—»2(U, 1) 


1.94(0.07) 


-9.2(0.1) 


3.3(0.2) 


0.55 


/""ITT rMT O/O O/O 

CH3UH 3(2,2)— »2(2, 2) 


1.6(0.4) 


-12(1) 


9(4) 


0.17 


OH3UH 3(1,3)— >2(1, 3) 


1.3(0.4) 


-10(1) 


8(3) 


0. 15 


O1T3 L^il (z,4 ) - — >4(z ,4 ) 


2.6(0.5) 


-9.4(0.5) 


5(1) 


0.49 


kj n;5Un 01 U, J, ) — >4^U, 1 J 


a . y ^ u . z j 


-8.8(0.1) 


4.1(0.6) 


U.DD 


v_ n 3 un □ ^ 1 ,4j — >4 ^ 1 ,4 j 


9 1 0^ 


-9.2(0.1) 


3.1(0.2) 


U.D4 


PHoDH 8(1 4~l >7(0 1\ 

V_J 1 1 3 11 O ^ 1 , 4: ) t i ^U,0 J 


n n(n 
£1 . \j y u . j 


-9.7(0.5) 


7(1) 


27 


CH3OH 11(1,4)— >10(2, 4) 


n 1 k(o oa) 


-9.9(0.8) 


6(2) 


0.02 


H 9 C O 3 1 2 — ^ ^11 


5.4(0.3) 


-9.4(0.1) 


2.2(0.1) 


2.2 




5.9(0.3) 


-9.8(0.1) 


5.5(0.2) 


1.0 


H2CO 2i2 — >ln 


8.1(0.6) 


-9.7(0.1) 


2.4(0.1) 


3.1 




3.1(0.6) 


-9.5(0.2) 


6(1) 


0.5 


H 2 13 CO 2i2 — In 


0.24(0.05) 


-9.4(0.4) 


4(1) 


0.05 


CH 3 CN 5(0)— 4(0) 


0.17(0.02) 


-10.2(0.2) 


2.6(0.4) 


0.06 


CH 3 CN 5(1)— 4(1) 


0.16(0.02) 


-9.9(0.2) 


2.9(0.4) 


0.05 


CH 3 CN 5(2)— 4(2) 


< 


0.07 K kms" 1 with Av=3.2 km s -1 




CH 3 CN 5(3)— 4(3) 


< 


0.07 K kms -1 with Av=3.2 km s -1 




CN 2—1 


6.0(0.3) 


-9.0(0.1) 


2.9(0.2) 


1.9 


N 2 D+ 3—2 


a 


=0.11 K kms -1 en 


Av= 1.5 km s -1 




DCO+ 2 — 1 


0.62(0.04) 


-9.58(0.04) 


1.5(0.1) 


0.38 


DCO+ 3—2 


0.76(0.07) 


-9.4(0.1) 


2.2(0.2) 


0.33 


D 2 CO 4q4- 3 03 


0.60(0.2) 


-9.1(0.5) 


4(1) 


0.15 



a In case of optically thin emission, r m cannot be determined and is set arbitrarily to 0.1. 



arising in a different region that the (1,1) and (2,2) lines. 
This interpretation is reinforced by the NH3 rotational di- 
agram (see Fig. 3). The three lines detected in NGC 7129- 
FIRS 2 cannot be fitted by one single straight line which 
implies the existence of at least two gas components with 
different rotation temperatures, a cold component traced 
by the (1,1) and (2,2) lines and a hot component only de- 
tected with the (3,3) line. Using the (1,1) and (2,2) lines, 
we derive a rotation temperature, Ti2 = 13 K and a column 
density, N(NH3)=4.9 10 14 cm -2 for the cold component. 
Because of the lack of radiative transitions between differ- 
ent K-ladders, the ammonia inversion lines are good ther- 
mometers of dense clouds. In fact, detailed radiative trans- 
fer calculations for N H 3 show that T fc =Ti 2 for T fc <20 K 
l|Danbvet al. l fl988). This low value of the kinetic tem- 
perature in the cold envelope implies that depletion could 
be important in this young protostar. 

Since we have not detected the (4,4) line, we can only 
derive lower and upper limits for the rotation temperature 
of the hot component. The lower limit is given for the exci- 
tation temperature between the (1,1) and (3,3) lines, and 
the upper limit is given for the excitation temperature be- 



tween the (3,3) line and the upper limit to the (4,4) line. 
We obtain that the hot component has a rotation temper- 
ature, T rot ^30 -90 K, which implies a lower limit of 50 
K to the kinetic temperature of this component. We have 
estimated N(NH3)^1.5 10 13 cm -3 for the hot component. 
This value has been derived from the integrated intensity 
of the (3,3) line assuming T rot =90 K and LTE conditions 
and is, in fact, a lower limit to the actual column density 
of the hot gas. 

The same procedure has been repeated for the am- 
monia lines towards LkHa 234. Similarly to NGC 7129 - 
FIRS 2, the ammonia emission towards LkHcv 234 cannot 
be fitted with one single rotation temperature (see Fig. 3). 
However, in this case the linewidth of the (3,3) line is more 
similar to those of the (1,1) and (2,2) lines (see Table 3). 
The lincwidths seem to increase monotonically with the 
energy of the transition, but there is not a jump between 
the linewidth of the (3,3) line and those of the others as 
in the case of the protostar. We have fitted a two compo- 
nents model to the rotational diagram in LkHcv 234 and 
obtained rotation temperatures of ~ 22 K and ~49-134 
K for cold and hot components respectively. This implies 
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7 





NGC 7129 - FIRS 2 








LkHa 234 






Molecule 


Comp 


i-rot (K.) 


IN (cm J 




Molecule 


Comp 


T rot (K) 


N 6 (cm" 2 ) 




NH 3 


Cold 


13 


4.9 10 14 




NH 3 


Cold 


22 


4.0 10 


8 




Warm 


31 - 87 


~ 1.5 10 13 






Warm 


49 - 134 


~ 1.1 10 13 




N 2 H+ 




13 a 


3.8 10 13 


-21" 


N 2 H+ 




22 a 


1.0 10 13 




H 13 CO+ 




13 a 


2.2 10 12 




H 13 CO+ 




22 a 


1.7 10 12 




CH3OH 




17 


7.8 10 14 




CH 3 OH p 




24 


1.4 10 14 






Hot core 


>80 


~ 2.0 10 14 








>250 


~ 4.0 10 14 














CH 3 OH e 




59 


4.8 10 14 




H 2 CO 


Narrow 


~ 10 


2.2 10 13 




H 2 CO 




~ 11 


8.0 10 13 






Wide 


~ 9 


4.8 10 13 














H 2 13 CO 




~ 10 a 


2.3 10 12 




H 2 13 CO 




~ ll a 


2.4 10 12 




CH 3 CN 


Hot core 


63 


3.6 10 12 




CH 3 CN 




<53 


> 2.0 10 12b 




CN 




5 


4.8 10 13 




CN 




6.5 


6.1 10 13 




HCN 




5 a 


1.6 10 13 




HCN 




6.5 a 


2.1 10 13 




N 2 D+ 




13 a 


5.4 10 11 




N 2 D+ 




22° 


< 2.3 10 11 




DCO+ 




8.5 


1.8 10 12 




DCO+ 




17 


4.2 10 11 




D 2 CO 




10 a 


2.0 10 12 




D 2 CO 






3.7 10 12 





a Assumed rotation temperature. 

p Assuming a point source. 

e Assuming a beam filling factor of 1. 



a kinetic temperatu re of Tfc~28 K for the cold envelope 

ijDanbv et al. Ill988|) and a lower limit of T fc >100 K for 

the hot component. Contrary to FIRS 2, the (1,1) line is 

moderately thick in this source. This allows us to estimate 

the size of the NH3 emission. Assuming that the exefta- 2 

tion temperature of the (1,1) line is equal to the rotation 1 

temperature of the cold component, we obtain a size of w 

8", which is lower than the size of the clump in the 1.3miu8 

continuum emission. 6 

0.4 



5.2. NoH+, H 13 CO + 



Recent studies in pre-stellar cores reveal that the widely 
used tracers of dense gas, C 18 and CS, are not adequateo 
to trace dense cold clumps. Both are strongly depleted- 1 
in the core at densities of a few 10 cm" " 3 and T fe <20^- 2 
llTafalla et, al. 112004) . Only nitrogenated molecules seem- 5 
to be unaffected by depletion. In fact, N2H+ seems to 
keep a constant abundance through the core and is an 
excellent tracer of the cold gas (e.g. ICaselli et al."lll999l 5 
iTafalla et al~ll2002l) . 

We have made small maps in the N2H + 1-0 and the 
H 13 CO+ 1-0 lines toward FIRS 2 and LkHa 234. In addi- 
tion we have observed the HC 18 + line towards the star 
position to have an estimate of the opacity of the H 13 CO + 
1-0 line. In Fig. 4 and 5 we show the N2H + 1-0 and the 
H 13 CO + 1-0 line integrated intensity maps together with 
the continuum map at 1.3mm, and in Table 2 and 3 we 
show the fits to the molecular lines. Due to the splitting 
of the N2H+ 1— >0 line, we can have an estimate of the line 
opacity using the same method as in the case of NH3. 

Intense centrally-peaked emission is observed in the 
continuum map at 1.3mm toward FIRS 2. The same mor- 
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Fig. 2. Spectra of the NH 3 (1,1), (2,2), (3,3) and (4,4) 
inversion lines towards NGC 7129 - FIRS 2 (left) and 
LkHa 234 (right). 

phology is observed in the H 13 CO + and N2H+ maps, al- 
though the profile of the H 13 CO + emission is flatter than 
that of N2H" 1 ". To quantify this difference in the emission 
profile, we have calculated the line intensity ratio, r\ — 
I m (N 2 H+)/I(H 13 CO+), where I m (N 2 H+) is the inten- 




Fig. 3. Rotational diagram of NH 3 in NGC 7129 - FIRS 2 and LkHa 234. The long-dashed line corresponds to the 
fit assuming a single rotation temperature for all the transitions. The solid and dotted lines are the warm and cold 
component respectively of the two-components model. 



sity of the main component of the N2H+ 1— >0 line and 
I{H 13 CO+) the intensity of the H 13 CO+ l->0 line, in a 
radial strip at 0" offset in declination. The ratio rj changes 
from 2.7 at the center of the clump to 1.5 at an offset of 
(-30",0"). Using the LTE approximation with T rot = 13 
K (based on our NH3 calculations) and assuming opti- 
cally thin emission in the H 13 CO + line, we derive that 
the N2H + /H 13 CO + abundance ratio changes from ~ 17 
in the (0,0) position to ~ 5 at (-30", ). Thus, the abun- 
dance of H 13 CO + relative to N2H + seems to decrease by 
a factor of ~ 3 towards the clump center. However, this 
change in the estimated H 13 CO + abundance could be due 
to the larger opacity of the H 13 CO + line. To constrain 
the opacity of the H 13 CO + 1— >0 line, we have observed 
the HC 18 0+ l->0 line towards the (0,0) position. The 
H 13 CO+ 1-vO /HC 18 0+ l->0 line intensity ratio is ~ 11, 
showing that the H 13 CO + 1-0 line is optically thin at this 
position. 

In the above calculations, we have assumed that the 
rotation temperature is uniform in all the strip. This as- 
sumption could be unrealistic since the density and kinetic 
temperature are expected to decrease with the distance 
from the star. We have used an LVG code to investigate 
if the different values of r± are due to a gradient in the 
physical conditions in the clump or/and are the result of a 
change in the relative abundance of both molecules. Since 
we are comparing the ground state lines of both species, 
the ratio r\ is very little dependent on the kinetic temper- 
ature and depends mainly on the hydrogen density and 
the N 2 H + /H 13 CO + abundance ratio. In Fig. 6 we show 
the ratio r\ and the opacity of the main hyperfine N2H" 1 " 
1— >0 line for a wide range of physical conditions assum- 
ing X= N 2 H+/H 13 CO+=3,7,15. The values of the opacity 
and r\ measured towards the star position can only be fit- 
ted assuming X=15 (see Fig. 6). A lower value of X would 
imply that the N2H + line is optically thin in contradiction 



with our observational results. On the contrary, the value 
measured at the offset (-30", 0") can only be fitted with 
X< 7. Thus we conclude that the gradient in the value of 
r\ cannot be due to the expected gradient in the excita- 
tion temperature of the observed lines across the clump, 
but to a gradient in the N2H + /H 13 CO + abundance ratio. 

As commented above, detailed studies in pre-stellar 
clumps show that the abundance of N2H 4 " remains con- 
stant in these cold clumps while H 13 CO + could suf- 
fer from depletion in the densest part l|Lee et al. I l2003t 
ICaselli et al . "1 12002^ . Assuming that N 2 H+ has a constant 
abundance in the clump, we need to assume an H 13 CO + 
depletion factor ^2 to fit our observations. 

The integrated intensity maps of the H 13 CO + 1— >0 and 
N2H + 1— >0 lines towards LkHa 234 are shown in Fig. 5. 
The N2H + 1— >0 map presents a very different morphology 
to that of the continuum and the H 13 CO + maps. While 
the continuum map presents a intense point-like emission 
at the position of LkHa 234, the N2H 4 " map peaks at 
the offset (-6", +18"). The map of the H 13 CO+ line show 
both peaks. As we will argue in the following, these differ- 
ent morphologies can only be explained if one assumes a 
different temperature and chemistry between both peaks. 

In Table 4 we show the estimated N 2 H+ and H 13 CO+ 
column densities in LKHa 234 assuming a rotation 
temperature of 22 K derived from the NH3 obser- 
vations. We have obtained N(N 2 H + )~1.0 10 13 cm" 2 
and a N 2 H+/H 13 CO+ abundance ratio of -6. The 
N 2 H+/H 13 CO + ratio in LkHa 234 is lower than that 
found in FIRS 2. 

We have repeated the same calculations for the off- 
set (-6", +18"). Continuum observations suggests that this 
clump is colder than that associated to LkHa 234. In fact, 
if we assume that the N2H" 1 " abundance is uniform in the 
region and a typical dust temperature of 30 K towards 
LkHa 234, we need to assume a dust temperature ~10 K 
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at the offset (-6", +18") in order to explain the measured 
continuum flux. Thus, we have assumed a rotation lower 
temperature, T rot — 10 K, in our LTE caculations. With 
these assumptions, we obtain a N 2 H + column density of 
-1.3 10 13 cm" 2 and a N 2 H+/H 13 CO+ abundance ratio of 
—10. This value is intermediate between those measured in 
LkHa 234 and FIRS 2. However, the difference is of a fac- 
tor of 2 which is within the uncertainty involved in these 
kind of calculations. Based on our results, we can propose 
an evolutionary trend based on the N2H + /H 13 CO + ratio. 
This ratio is maximum in the IM Class object FIRS 2 
where molecular depletion is significant, it may take an in- 
termediate value in the infrared low-mass star IRS 6 and 
is minimum in the envelope of the more evolved object, 
the HBe star LkHa 234. 

5.3. CO outflows, SiO 

The SiO abundance is strongly enhanced (up to several or- 
ders of magnitude) in shocked regions. How ever, its abun- 
dance is very low in dark clouds and PDRs l|Schilke et al. I 
l200lh . Because of this peculiarity, SiO is used as a diag- 
nostic for shocks in both galacti c and extragalactic re- 
gions (e.g. iBachiller et al. j jl99lt I Martm- Pintado et al. I 
ll992HOarcia-Burillo et al. 11200(1 EoOll) . In particular, it 
is commonly used to look for the energetic outflows asso- 
ciated with the youngest stellar objects. We have carried 
out maps of the J=2^1 and 3^2 rotational lines of SiO 
around FIRS 2 and LkHa 234 to study the physical condi- 
tions and chemistry of the bipolar outflows found in these 
regions. 

Interferometric and single-dish CO 2— >1 observations 
show the existence of two bipolar outflows associated with 
FIRS 2. The axis of these outflows form an angle of almost 
90° giving a quadrupolar morphology to the spatial distri- 
bution of the high velocity gas of the region. In Fig. 7 we 
show the m aps of the high veloc ity 12 CO 2— >1 emission 
reported bv iFuente et al. I l)200l[) . Interferometric obser- 
vations shows that the NE-SW outflow, hereafter FIRS 
2-out 1, is the one associated with the millimeter source 
FIRS 2-MM1. An intense blue lobe is detected in the CO 
2— >1 emission in this outflow while the red lobe is only 
marginally detected. The NW-SE outflow, hereafter FIRS 
2-out 2, seems to be assoc iated with an infrared star FIRS 
2-IR l|Fuente et al. Il200ll) . In FIRS 2-out 2, the red lobe 
is more prominent than the blue one, giving a peculiar 
appearance to the high velocity gas around FIRS 2. The 
peaks of the 12 CO emission in this lobe present a "Z"-like 
shape which suggests that the axis of the outflow FIRS 
2-out 2 is precessing. Precession has been fo und in bipolar 
outflo ws associated to young low-mass stars l)Cueth et al. I 

US- 

In Fig. 8 we present the integrated intensity map of the 
SiO 2— >1 line towards FIRS 2. The SiO emission is only ob- 
served along the axis of the two outflows. Furthermore the 
SiO peaks are not located at the star position but in the 
bow shocks located at the end of the outflow lobes. Thus, 



the SiO emission seems to be closely related to the bipo- 
lar outflows. The spectra of the SiO lines towards some 
selected positions are shown in Fig. 8. The SiO 2— >1 line 
towards the star position presents two velocity compo- 
nents, a narrow component with Av — 1.6 km s -1 which 
is centered at the velocity of the ambient cloud, V/ sr — 

— 10.0km s _1 , and a wide component, Aw— 7.0 km s^ 1 
which is centered at V/ sr - — 7.0+1.0 km s" 1 . The compo- 
nent at —7.0 kms -1 corresponds to a well defined high 
velocity clump which peaks at the position (-10", -10") 
(see Fig. 8). Hereafter we will refer to this clump as Rl. 
There exists a counterpart blue clump at a velocity, V; sr - 

— 12.5+0.5 km s" 1 which peaks at the position (50", 50") 
(hereafter Bl). At this position, the profile of the SiO 2— >1 
line only presents a wide component with Aw— 7 km s _1 . 
The high velocity clumps Rl and Bl have well defined ve- 
locities a nd positions like the "bullets" found in low-mass 
stars (see IBachiller lfl996|) . The jet- like morphology of the 
SiO emission along the outflow axis as well as the existence 
of bullets argues in favor of the youth of this outflow. 

The component at the velocity of the ambient cloud, 

Vi sr 10.0 km s -1 , is also present in FIRS 2-out 2. The 

emission of this narrow component in FIRS 2-out 2 is sur- 
rounding the red lobe as detected in 12 CO. In fact, the 
narrow component is located adjacent to the peaks of the 
high velocity CO emission suggesting that SiO emission 
is tracing the molecular cloud gas entrained in the out- 
flow. Narrow SiO components widespread in the molec- 
ular cloud in which the bipolar outflow is embedde d has 
been detected in other Class protostars l|Codella et, al. I 
1999). We propose an interpretation in which the mor- 
phology of the SiO emission is related to the evolutionary 
stage of the outflow. In FIRS 2-out 1, the SiO emission 
has a jet-like morphology and is concentrated in "bullets" 
ejected by the exciting star. This "jet-like" morphology is 
als o observed in the in terferometric 12 CO image reported 
by l|Fuente et al. l200lh . The SiO emission in FIRS 2-out 2 
is surrounding the red CO lobe. We propose that in this 
case the SiO emission traces the material adjacent to the 
cavity walls excavated by the outflow which is being en- 
trained into the outflow. The different profiles of the SiO 
emission in FIRS 2-out 1 and FIRS 2-out 2 are clearly 
seen in Fig. 8. 

We have calculated the SiO 3— >2/2^1 line intensity 
ratio, r32, at the (0",0") position by degrading the an- 
gular resolution of the SiO 3^2 map to that of the SiO 
2— >1 one. A value of r32 — 0.8 is found in the three chan- 
nels centered at the ambient velocity, while a value 
> 1.0 is found for the high velocity gas. This reveals a 
higher excitation temperature for the high velocity gas. 
In order to determine the physical conditions for both 
components, we have used an LVG code. Since sug- 
gests different physical conditions for the different velocity 
ranges, we have carried out LVG calculations for the ambi- 
ent, moderate velocity and high velocity gas components 
(see Table 6). Taking into account the kinetic tempera- 
ture derived from the NH3 lines we have assumed Xfc=15 
K for the ambient narrow component and Tfe=50 K for 
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Table 5. Beam averaged N2H + and H 13 CO + column densities 



Source 


Trot (K) 

(K) 


N(H"CO+) 

(cm- 2 ) 


N(N 2 H+) 
(cm" 2 ) 


N 2 H+/H i3 CO + 


FIRS 2 


13 


2.2 10 12 


3.8 10 la 


17 


(-30,0) 


13 


8.0 10 11 


4.2 10 12 


5 


LkHa 234 


22 


1.7 10 12 


1.0 10 13 


6 


(-6.+18) 


10 


8.8 10 11 


8.5 10 12 


~ 10 




Fig. 4. Continuum map at 1.3mm and integrated intensity maps of the observed lines in NGC 7129 - FIRS 2. The 
box size is 100" x 100" and is centered at the star position. Contour levels are: a) 39.2 mJy/beam, and from 78.5 to 
706.2 mJy/beam by 78.5 mJy /beam ;b) 0.19 to 1.7 by 0.19 K kms" 1 ; c) 0.2 to 8.2 by 0.8 K kms" 1 . d) 0.25 to 4 by 
0.5 K kms- 1 ; c) 1.5 to 4 by 0.5 K kms" 1 ; f) 1.5 to 13.5 by 1.0 K kms" 1 ; g) 2 to 20 by 2 K kms" 1 ; h) 0.1 to 0.5 by 
0.1 K kms" 1 ; i) 1 to 14 by 1 K kms" 1 ; j) 1 to 21 by 2 K kms" 1 ; k) 2 to 7 by 1 K kms" 1 ; 1) 0.2 to 1.9 by 0.2 K kms" 1 ; 
m) 0.8 to 3.4 by 0.4 K kms" 1 ; n) 0.2 to 3 by 0.4 K kms" 1 . 
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Fig. 5. Same as Fig. 1 for LkHev 234. Contour levels are: a) 39.2 mJy/beam, 78.5 to 706.2 mJy/beam by steps of 78.5 
mJy /beam ; b) 0.19 to 1.7 by 0.19 K kms" 1 ; c) 0.2 to 8.2 by 0.4 K kms" 1 ; d) 1.5 to 15 by 1.5 K kms" 1 ; e) 1 to 11 
by 1 K kms" 1 ; f) 1.5 to 16.5 by 1.5 K kms" 1 ; g) 2 to 11 by 1 K kms" 1 ; h) 0.5 to 5 by 0.5 K kms" 1 ; i) 1 to 6 by 1 K 
kms -1 



the high velocity components. The derived densities and 
column densities are shown in Table 6. We are aware that 
the NH3 and SiO emissions could arise in different regions 
and the assumed kinetic temperatures are very uncertain. 
But since the SiO lines are optically thin, the derived col- 
umn densities are weakly dependent on the assumed ki- 
netic temperature and are accurate within a factor of ^2. 
The values of Table 6 clearly show that the SiO abundance 
is larger by almost 2 orders of magnitude in the high ve- 
locity gas of FIRS 2-out 1 than in the narrow component 
associated with FIRS 2-out 2. Within FIRS 2-out 1, we 
also detect a gradient in the SiO abundance, being the 
SiO/ 13 CO ratio 2 orders of magnitude larger in the high 
velocity component than in the ambient component. This 
enhancement of the SiO abundance in the high velocity 
gas is als o observed in very young: lo w-mass protostellar 
outflows ( Bachill er et al. I ll991. 20 0l[) and interpreted as 
due to the chemical gas processing by the energetic shocks 
associated with the high velocity "bullets" . 

Hydrogen densities are also quite independent on the 
assumed temperature for Tfc>50 K. The estimated hydro- 
gen density decreases by only a factor of ^4 if we change 
the kinetic temperature from 15 K to 100 K. Thus, the 
estimated hydrogen densities are accurate within a factor 
of ~4. Our calculations show that the density seems to in- 
crease from a few 10 5 cm" 3 in the ambient component to 



>10 6 cm -3 in the high velocity gas at the star position. We 
have also carried out LVG calculations for Rl and Bl, and 
at the offset (50", -50") where only the narrow component 
has been detected. The density in the bullets Rl and Bl 
is <Jl0 5 cm -3 while the density in the narrow component 
detected at the offset (50",-50") is of a few 10 4 cm" 3 . This 
density is an upper limit to the hydrogen density in the 
narrow component, since the kinetic temperature is never 
expected to be lower than 15 K. Thus, we consider that the 
estimated difference in the hydrogen density of the wide 
and narrow components is reliable. The lower density of 
the narrow component also supports our interpretation of 
the SiO emission arising in the gas of the molecular cloud 
surrounding the outflow. 

We have observed a map of 80"x80" around LkHa 
234 in the SiO 2-^1 and 3^2 lines. We have detected SiO 
emission only towards the star position. Since we have in- 
tegrated double time towards the star position than in the 
other map positions, we cannot exclude the possibility of 
SiO emission at the same level in other positions. Thus 
we have poor information about the spatial distribution 
of the emission. Regarding the kinematical information, 
the large linewidth of the SiO lines, Av>3 km s _1 , com- 
pared to that of the (1,1) and (2,2) ammonia lines sug- 
gests that the emission arises in the warm component. 
The weakness of the SiO 2— >1 line emission as well as the 
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Fig. 6. Plots of the r% = I(N2H + )/ I(H 13 CO + ) ratio and opacity of the main N2H + hyperfinc component (solid 
contours), and r m (dashed contours), as a function of the H 13 CO+ column density, N(H 13 CO + ), and hydrogen 
density, n, for a kinetic temperature of Tk = 15K and an N2H + /H 13 CO + abundance ratio, X=3,7,15. We have shaded 
with dark grey the region with n = 1.5±0.3 ( ratio measured at the 32" offset from the star position) and with a 
clear color r\ = 2.7±0.6 (ratio measured at the star position). We need a depletion of at least a factor of 2 towards 
the center in order to fit both tau m and n in both positions. Contour levels are: a) r± = 0.8, 1, 1.2 and r m = 0.1, 0.5, 
1.0, 2, 5,10; b) n = 1.8, 2.2, 2.7, 3.2 and r m = 0.1, 0.5, 1.0, 2, 5, 10; c) n = 3.2, 5, 7 and r m = 0.1, 0.5, 1.0, 2, 5,10 



lack of information about the source size make any es- 
timate of the hydrogen density very uncertain. We have 
derived the SiO column density assuming n^ 5 10 5 cm" 3 . 
With this assumption we derive a SiO column density of 
^lO 11 cm~ 2 . This value is lower by a factor of >6 than the 
total SiO column density towards FIRS 2. The SiO/ 13 CO 
ratio towards this st ar is similar to that found in PDRs 
iScnilke et a,l.ll200l. 



Table 7. SiO emission 



In Table 7 we present a summary of the SiO observa- 
tions. There is a clear evolutionary trend in the SiO behav- 
ior. The youngest outflow FIRS 2-out 1 presents intense 
SiO emission at high velocity with SiO abundances as high 
as ^10~ 8 . Towards the more evolved outflow FIRS 2-out 2 
we have detected only a weak SiO line at ambient velocity. 
The SiO abundance in this component is ^10~ 10 , i.e., two 
orders of magnitude lower than the SiO abundance in the 
high velocity gas associated with FIRS 2-out 1 but still 
larger than the SiO abundance in PDRs and dark clouds. 
Towards LkHa 234 we have detected SiO emission at the 
ambient velocities with a fractional abundance of 10~ 12 . 
This abundance is similar to that measured in PDRs and 
could be associated with the PDR produced by LkHa 234. 

This evolutionary trend confirms that SiO is a good 
tracer of energetic shocks. The SiO abundance is highly 
enhanced when the shocks are strong enough to re- 
lease the silicon from the grains to the gas phase 
I Martin-Pintado et al. Il992j) . In a MHD shock model the 
rele ase of Si to the gas phase requires V K hnr.k> 4 km 
s -1 ((Flower. Pineau des Forets. Field, fc Mnvll !)!)(!). This 
is consistent with the trend of having larger SiO abun- 
dances in the higher velocity gas. Because of projection 
effects, the velocity we measure is a lower limit to V shock- 
As the protostar evolves, the ouflow fades and the amount 
of high velocity molecular gas decreases. This produces a 



Outflow 


Age (Myr) 


Component 


X(SiO) 3 


FIRS 2-out 1 


0.003 1 


High velocity 


~ 10" 8 


FIRS 2-out 2 


0.005 1 


Ambient 


~ icr 10 


LkHa 234 


~ 0.1 2 


Ambient 


~ 10" 12 -10 -11 



Estimated age of the outflow assuming a velocity of 15 km 
s _1 and the lobe sizes measured from Fig. 1. 
2 Stellar age from lFuente et al. I |[l998) 

Fractional abundance estimates assuming X( 13 CO) = 2 
1(T 6 and X(H 13 CO+)=l(T 10 



decrease in the SiO abundance. When the bipolar molec- 
ular outflow stops, the SiO abundance around the star 
decreases to the typical value in PDRs. 

5.4. CS, C 34 S 

We have carried out a map of the CS 3—^2 emission to- 
wards FIRS 2 and LkHa 234. In addition we have observed 
the C 34 S 2->l, 3^2 and 5^4 lines towards the (0,0) posi- 
tion. The spectra towards the (0,0) position are shown in 
Fig. 9 and the Gaussian fits are shown in Tables 2 and 3. 
The integrated intensity maps of the CS J=3^2 line are 
shown in Fig. 4. 

Similarly to the case of SiO, the profiles of the CS and 
C 34 S lines towards FIRS 2 present high velocity wings 
at red-shifted velocities. However, the terminal velocity 
of these wings is lower than those of the SiO lines. The 
velocity range of the CS emission is between -14 and -5 
km s _1 , i.e. , it is only detected at the ambient velocities 
and in the moderate velocity component of the SiO emis- 
sion. The spatial distribution of this component in the CS 
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100 50 -50 -100 

Aa(arcseconds) 



Fig. 7. High velocity CO emission in FIRS 2. The blue-shifted velocities are shown in the top panels, and the red- 
shifted velocities in the bottom panels. Straight lines indicate the outflow axis. In the case of FIRS 2-out 2 we have 
drawn two axis because of outflow precession. Contour levels are: a) 1 to 6 by 1 K kms -1 ; b) 2 to 22 by 2 K kms -1 ; 
c) lev 7 to 70 by 7 K kms" 1 ; d) 7 to 80 by 7 K kms" 1 ; c) 4 to 50 by 4 K kms" 1 ; f) 1 to 8 by 1 K kms" 1 . 



3^2 line is different from that of the same component of 
the SiO emission, suggesting that both molecules might 
be tracing different gas components even when comparing 
the same velocity range. 

We have carried out LVG calculations to estimate the 
physical conditions of the gas emitting in CS. Different line 
ratios are found at the ambient velocities and in the mod- 
erate velocity range in the C 34 S lines (see Fig. 9). Thus, we 
have estimated the physical conditions in the two velocity 
intervals. We have derived a density >3 10 6 cm -3 for the 
ambient component in FIRS 2. This density is larger by a 
factor ~6 than the density estimated from SiO lines. This 
difference is not due to the assumed kinetic temperature, 
since even assuming a kinetic temperature of 100 K the 
density would decrease only by a factor of ~4. However, 
the densities derived from the C 34 S lines in the moderate 
range are similar to those derived from the SiO lines. We 
propose that at least part of the C 34 S lines emission at 
ambient velocities does not arise in the outflow. We spec- 
ulate with the possibility of the C 34 S emission arising in 
a hot core. 

The CS and C 34 S lines observed in LkHa 234 are 
shown in Fig. 9 The line-widths of the C 34 S lines are sim- 
ilar to that of the SiO and the NH 3 (3,3) lines suggest- 
ing that they are tracing a warm component. The emis- 
sion of the CS lines is concentrated towards the star po- 
sition, therefore the connection between the outflow and 
the CS emission is not clear. We have derived a density 
~10 5 cm~ 3 for this warm component. This density is lower 
than the typical density of the hot cores associated with 



massive stars. Since LkHa 234 is a visible star, the enve- 
lope has already been disrupted by the star and the UV 
radiation is escaping through the envelope. A dense (n> 
10 6 cm -3 ) hot region similar to the hot cores in massive 
stars is not expected at this evolutionary stage. 

5.5. CH 3 OH 

We have observed up to 17 methanol lines towards the 
(0,0) position in FIRS 2 and LkHa 234. The large number 
of lines allows us to have a good estimate of the rotation 
temperature and the CH3OH column density. In addition, 
we have obtained a small map in the intense CH3OH 2— >1 
and 5^4 lines towards FIRS 2. The CH3OH maps around 
FIRS 2 are shown in Fig. 4. The CH3OH emission peaks 
towards the bullet Rl suggesting that these lines arise 
mainly in the molecular outflow. 

The high- velocity resolution spectra of the CH3OH 
lines are shown in Fig. 10. The velocity profiles of the 
CH3OH lines towards FIRS 2 are similar to those found 
in the SiO lines. The low energy lines present two velocity 
components, a narrow component at the ambient velocity 
and a wide wing which extends to red-shifted velocities 
and is centered at -6±1 km s _1 . In higher energy tran- 
sitions, the narrow component becomes weaker and only 
the wide red-shifted one is detected. This is consistent 
with the integrated intensity maps of the CH3OH 2— >1 
and 5^4 lines shown in Fig. 4 which show a strong peak 
emission at the position of the bullet Rl. In order to fur- 
ther investigate the nature of the CH 3 OH emission, we 
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Table 6. LVG calculations : SiO 



Source 


Vel. 
(km 


range 
s" 1 ) 


lfc(K) 
(K) 


n 

(cm' 3 ) 


IN ( LU ) 

(cm- 2 ) 


TVT fQli~\\ 

IN (blU ) 

(cm- 2 ) 


IN [ti ) 

(cm -2 ) 


SiO/ 13 CO 


SiO/H 13 CO+ 


FIRS 2 (Star) 


-13.5 


- -11.5 


50 


5 10 b 




1.5 10 11 


1.0 10 11 




1.5 




-11.5 


- -8.5 


15 


5 10 5 




1.1 10 12 


1.5 10 12 




0.7 




-8.5 - 


-4.5 


50 


4 10 5 




1.2 10 12 










-4.5 - 


-0.5 


50 


1 10 6 




4.7 10 11 










Total 










2.8 10 


1.5 10 12 




1.9 


(-10" ,-10") Rl 


-13.5 


- -11.5 


50 


5 10 5 * 


3 10 15 


4 10 11 




0.0001 




(FIRS 2-out 1) 


-11.5 


- -8.5 


15 


4 10 5 


4 10 16 


3 10 12 


1 10 12 


0.00007 


3 




-8.5 - 


-4.5 


50 


4 10 5 


3 10 15 


1.5 10 12 




0.0005 






-4.5 - 


-0.5 


50 


3 10 4 


o -I nl4 

3 10 


1.5 10 12 




0.005 






Total 








4.6 10 16 


6.4 10 12 


1 10 12 


0.0001 


6 


(50",50") Bl 


-17.5 


- -12.5 


50 


2 10 5 


1.5 10 15 


1.5 10 12 




0.001 




(FIRS 2-out 1) 


-12.5 


- -8.5 


50 


2 10 5 


5.5 10 16 


1.5 10 12 




0.00003 






-8.5 - 


-4.5 


50 


2 10 5 


2.0 10 15 


2.0 10 11 




0.0001 






Total 








5.8 10 16 


3.2 10 12 


< 2 10 12 


0.00005 


> 1 


(50,-50) 02 


-11.5 


- -8.5 


15 


8 10 4 




1.0 10 11 


3 10 11 




3 


(FIRS 2-out 2) 




















LkHa 234 
(IRS 6-out 1) 


-12.5 


- -8.5 


50 


5 10 5 * 


1.3 10 17 


2.0 10 11 


1.0 10 12 


0.000001 


0.2 



* Assumed value for the hydrogen density 



Table 8. LVG calculations : CS, C 34 S 



Source Vel. range T fc (K) n (cm' 3 ) N s (C 34 S)(c m - 2 ) Q s 



NGC 7129 - FIRS 2 


-10.5 


- -9.5 km s- 1 


15 


3 10 b 


8.0 10 11 


~9" 




-9.5 - 


-6.5 km s -1 


50 


6 10 5 


7.1 10 11 


^4" 




Total 








1.5 10 12 




LkHa 234 


-11.5 


- -9.5 km s- 1 


50 


2 10 5 


9.6 10 11 


14" 




-9.5 - 


-5.5 km s -1 


50 


4 10 5 


7.8 10 11 


6" 




Total 








1.7 10 12 





have fit the lines with Gaussian profiles and make some 
correlations. In Fig. 11 we plot the central velocity and 
lincwidths of the CH 3 OH lines vs. the upper state en- 
ergy of the observed transition. It is clearly seen that the 
lincwidths of the CH3OH lines increase with the upper 
state energy of the transition at moderate energies. In 
fact, they increase from ^3 km s _1 in the low energy 
transitions to ^7 km s _1 in transitions with E u >50 K. 
But this trend is not valid for higher energy lines which 
seem to have a constant linewidth of ~4 km s _1 . A similar 
behavior is found when one compares the velocity of the 
line with the upper state energy. The line velocity changes 
from v^ — 10 km s' 1 to ~— 7.5 km s' 1 when the energy 
increases from ~10 K to 50 K. However, for higher ener- 
gies, the line velocities seem to make the reverse way and 
change from ^—8 km s _1 to —10 km s _1 . This suggests 
that the emission of the low and moderate energy tran- 
sitions ( E u >50 K) arise in the molecular outflow. The 
correlation found between the linewidth and the energy of 
the transition suggests that the high velocity gas is asso- 



ciated with higher excitation temperatures. This result is 
consistent with the density estimated from the SiO lines. 
For E„> 50 K, there is a jump in the line velocity which 
go back to ~ -10 km s -1 . We propose that this could be 
due to the existence of a hot component in the CH3OH 
emission in addition to that related to the bipolar outflow. 

In Fig. 11 we show the CH3OH rotational diagram for 
FIRS 2. All the observed CH3OH transitions cannot be fit- 
ted with a single straight line. We need to assume at least 
two rotation temperatures to fit the observational data. 
The low energy transitions (E u <50 K) are well fit with 
a T rot ~17 K, while the high energy transitions require a 
higher rotation temperature, T rot >80 K. We propose the 
existence of a "hot core" component which dominates the 
emission in the high energy transitions. 

Several CH 3 OH lines have also been observed in LkHa 
234. In this case all the lines are centered at the ambi- 
ent velocity. However, there are important variations in 
the linewidths of the observed lines. Like in the case of 
FIRS 2, the linewidths seem to increase with the energy 
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of the upper level for E M <50 K (see Fig. 12). We have 
only detected the methanol lines towards the (0,0) posi- 
tion and consequently, we have no information about the 
size of the emitting region. Thus, we have considered the 
two limiting cases of a point source and a beam filling fac- 
tor of 1 to make the rotational diagram. In the first case, 
we need two gas components to fit all the observed tran- 
sitions. The cold one would have T rot ^24 K, and the hot 
one, T ro t>250 K. But if we assume a beam filling factor of 
1, all the oberved transitions are well fit with a T ro t~60 
K. With our data, we cannot discern between these two 
cases. 

5.6. H 2 CO 

We have observed two H2CO rotational transitions and 
one of the rarer isotope H 2 13 CO + toward FIRS 2 and 
LkHa 234. The obtained spectra are shown in Fig. 10. 
Similarly to other molecules, the H2CO spectra towards 
FIRS 2 present two well differentiated components, a nar- 
row one centered at ^—9.6 km s _1 with a linewidth 
of ~1.3 kms s _1 and a much wider one centered at 
^—9.0 km s _1 . However, the wide component has not the 
typical Rl profile observed in the SiO, CS and methanol 
lines at the star position. While the profiles of the SiO,CS 
and methanol lines present only red wings, the H 2 CO lines 
present a quite symmetric profile with blue and red wings. 
Consequently, the central velocity of the wide H2CO com- 
ponent is similar to that of the ambient gas and the 
linewidth is as large as Av=7 km s _1 . To further inves- 
tigate the nature of these components we have studied 
the integrated intensity maps of the H 2 CO 2i 2 — >ln line 
for the different velocity intervals (see Fig. 13). The most 
striking feature could be the jet-like morphology observed 
in the H 2 CO emission at blue velocities (from -15 to - 
10 km s _1 ). At red velocities, the emission is maximum 
at the offset (-7", -7") which is located close to the bullet 
Rl. Since the wide component have a very well differ- 
entiated profile, we have been able to subtract the wide 
component to the observed spectra and mapped both com- 
ponents separately. Our results are quite suggestive. The 
wide component presents a jet-like morphology with the 
maximum towards the position Rl. The morphology of the 
narrow component is an intense ridge which surrounds the 
jet. This strongly suggests that the narrow component is 
tracing the shocked gas of the molecular cloud which is in- 
teracting with the jet. The maximum of this narrow com- 
ponent coincides with the position where the bullet Rl 
impinges in the cloud. We have derived rotation temper- 
atures and column densities in the narrow and wide com- 
ponents separately, and obtained similar excitation condi- 
tions in both components. Thus, although the kinematics 
is clearly different, the physical condition of both compo- 
nents are quite similar. 

In Fig. 10 we show the H 2 CO spectra towards LkHa 
234. The profiles of the H 2 CO lines in this source, also 
suggest the existence of a narrow and wide components. 



However, these two component cannot be easily separated. 
For this reason, we have derived rotation temperatures by 
considering the sum of the two components. We obtain a 
rotation temperature and H 2 CO column density similar 
to those obtained in FIRS 2. In Fig. 5 we show the inte- 
grated line intensity maps in this source. Similarly to the 
case of FIRS 2 we find emission along the outflow and in 
a direction perpendicular to it. Thus far, no bipolar out- 
flow has been detected in this direction. Thus, this H 2 CO 
emission is associated with the flattened clump in which 
the Hcrbig Be star LkHa 234 is embebded that is being 
heated by the recently born star. 

5.7. CH 3 CN 

We have observed the CH 3 CN 5^4 and 13^12 lines to- 
wards FIRS 2 and LkHa 234. Because of the rotational 
structure of CH3CN, one can observe several lines at dif- 
ferent energies very close in frequency. This allows us to 
estimate the rotation temperature avoiding observational 
errors and the uncertainty due to the unknown source size. 
We have carried out these calculations towards our two 
sources. Unlike the other molecules observed, we do not 
detect a cold component in the CH3CN lines, but only 
the warm one. The detection of a hot CH3CN component 
with Tfc > 63 K in FIRS 2 shows the existence of a hot 
core in this object. 

CH3CN seems to be the best tracer of hot cores in 
these intermediate-mass stars. Contrary to CH 3 OH and 
H 2 CO whose low energy lines arise mainly in the bipolar 
outflow, the rotational lines of CH3CN seems to arise in 
the hot core and provide a good measure of the kinetic 
temperature of this hot component. 

5.8. CN,HCN 

We have mapped a small region around FIRS 2 and LkHa 
234 in the HCN l->0 and CN 1^0 and 2->l lines. The in- 
tegrated intensity maps around FIRS 2 are shown in Fig. 4 
and 5. Thus far, we have two different kind of molecules 
depending on the morphology of their emission in FIRS 2. 
The first group is formed by the molecular ions N 2 H + and 
H 13 CO + which peak at the YSO position and trace the 
bulk of the cold envelope. The second group is formed by 
CS, SiO, CH 3 OH and H 2 CO peaks at the position Rl and 
its emission has an important contribution from the molec- 
ular outflow FIRS 2-out 1. The morphology of the CN l->0 
emission is different from those of these two groups. The 
CN emission does not peak at the star position nor at the 
position of the bullet Rl, but to the north, suggesting a 
new chemical differentiation in the protostellar envelope. 
The morphology of the HCN emission is also different from 
those of the other molecular lines. In fact, the HCN emis- 
sion peaks to the NW from the star position, in a position 
intermediate between the bullet Rl and the peak of the 
CN emission. The large linewidths of the HCN 1— >0 line 
(see Fig. 15) show that, in contrast to the CN 1— >0 line, 
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Fig. 13. Panels a-d are the integrated intensity maps of the H2CO 2i2— ►In line towards NGC 7129 - FIRS 2 for 
different velocity intervals. In panel e) we show the integrated intensity map of the H2CO lines after subtracting the 
Gaussian fit to the wide component from all the spectra. In panel f) we show the area of this wide component as 
derived from the Gaussian fit. The axis of the outflow FIRS 2-out 1 and FIRS 2-out 2 are indicated. Contour levels 
are: a) 0.5 to 5 by 0.5 K kms" 1 ; b)1.0 to 6.5 by 0.5 K kms' 1 ; c) 0.5 to 5 by 0.5 K kms" 1 ; d)l to 14 by 1 K kms" 1 ; e) 
0.5 to 5 by 0.5 K kms" 1 ; f) 1 to 11.5 by 1 K km s" 1 . 



the emission of the HCN 1— >0 line has a significant con- 
tribution from the molecular outflow giving rise to this 
peculiar morphology. 

The radicals CN and HCN are known to be spe- 
cially abundant in PDRs. In particular, the CN/HCN ra- 
tio has been successfully used PDR tracer in dif- 
ferent kind of objects. In Fig. 16, we show the maps of 
the (CN l->0)/(HCN l->0) intensity ratio in FIRS 2 and 
LkHa 234. The (CN 1 -►())/ (HCN l->0) line intensity ra- 
tio is maximum at the star position and to the north, 
forming a conical feature with the star at its apex. We 
have been estimated the CN rotation temperature from 
the (CN 2->l)/(CN l->0) line intensity ratio (see Table 4). 
Assuming the LTE approximation and the same rotation 
temperature for CN and HCN, we obtain a CN/HCN 
abundance ratio of ^3 at the star position. This value 
is similar to those found in PDRs and suggests that the 
gas chemistry in this conical feature is being affected by 
the UV radiation from the protostar. However, the axis 
of this conical feature seems to be more similar to that 
of the outflow FIRS 2-out 2 than to that of the outflow 
FIRS 2-out 1. This suggests that the PDR traced by the 
high CN/HCN ratio could be related to the star driving 
the outflow FIRS 2-out 2 instead of to the Class object. 
In fact, the PDR could be formed in the walls of the cav- 
ity excavated by the outflow FIRS 2-out 2 when they are 
illuminated by the exciting star. But observations with 
higher spatial resolution are required to conclude about 
this point. 

We have also observed the CN and HCN lines towards 
LkHa 234. In this case, the linewidths of the CN and HCN 



lines arc similar, and in agreement with those found in the 
warm component. But, contrary to most of the observed 
molecular species, the CN and HCN emission do not peak 
at the star position but to the north forming a conical fea- 
ture. We have calculated the CN/HCN integrated inten- 
sity ratio in the region. Surprisingly, the CN/HCN ratio is 
minimum at the star position and maximum at the border 
of the clump as traced by the 1.3mm observations suggest- 
ing that the clump is illuminated from outside. Making 
column densities estimates, we derive N(CN)/N(HCN)~3 
at the star position. 

Thus, the CN/HCN fractional abundance ratio at the 
star position is equal (within the uncertainties) in FIRS 2 
and LkHa 234 and consistent with the expected value in a 
PDR. However the behavior of the CN/HCN ratio is very 
different in the rest of the envelope. In the case of FIRS 2, 
the CN/HCN ratio decreases outwards from the star as 
expected from a PDR illuminated from the interior and 
with an optically thick envelope. In the case of LkHa 234, 
the CN/HCN ratio is quite constant inside the clump and 
increases at the edges. This suggests that the clumps is 
also illuminated from outside (the clump is located at the 
border of an HII region). Since the envelope is less massive 
than that associated with FIRS 2, the whole envelope can 
be considered as a PDR. 

6. Deuterated compounds 

In order to derive the deuterium fractionation we have 
observed the DCO+ 2->l and 3^2, N 2 D+ 3^2 , and 
D 2 CO 4 4^3q3 lines toward the studied regions. In Fig. 15 
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we show the spectra of these lines toward FIRS 2 and 
LkHa 234. The integrated line intensity maps of the 
DCO+ 2-»-l line toward FIRS 2 are shown in Fig. 4. 

The linewidths of the DCO+ and N 2 D+ lines are -1.0 
km s _1 suggesting that they arise in the cold compo- 
nent of the envelope like the non-deuterared compounds 
HCO+ and N 2 H+. The linewidth of the D 2 CO line is -4 
km s _1 like those of the lines arising in the warm com- 
ponent, and in particular, the lines of the chemically re- 
lated species H2CO and H2 13 CO. In Table 9 we show the 
DCO+/H 13 CO+, N 2 D+/NH 2 + and D 2 CO/H 2 CO abun- 
dance ratios in both sources. The DCO + /H 13 CO + abun- 
dance ratio is a factor of 20 lower in LkHa 234 than in 
FIRS 2. This factor is so large that cannot be due to the 
H 13 CO + depletion but to a different value of the deu- 
terium fractionation in these cold envelopes. Thus, we 
propose that the deuterium fractionation in the cold en- 
velope decrease during the protostellar evolution. As we 
will comment in detail in next section, this increase in the 
deuterium fractionation can be understood as the conse- 
quence of the envelope warming during the protostellar 
evolution. 

A very different case is the D 2 CO/H 2 CO abundance 
ratio which increases by a factor —1.5 from FIRS 2 to 
LkHa 234. Since a factor of 1.5 is within the uncertaini- 
ties of our column density estimates, we conclude that the 
deuterium fractionation, as measured by the D 2 CO/H 2 CO 
abundance ration, seems to be constant (or slightly in- 
crease) in the warm component during the protostellar 
evolution. Thus, the evolution of the deuterium fraction- 
ation in the warm envelope seems to follow a different 
trend than in the cold envelope. The evaporation of the 
icy grain mantles is very likely the main responsible of this 
behavior. 

7. Discussion 

7.1. Physical structure of the YSOs envelopes and its 
evolution 

Our data show the existence of at least two well differenti- 
ated components in the envelope of FIRS 2 and LkHa 234, 
the cold envelope traced by the low energy lines of NH3 , 
N 2 H + and H 13 CO + , and the warm envelope traced by the 
CS, CH3OH and H 2 CO lines. These two components can 
be differentiated observationally by their kinematics, the 
morphology of their emission and by their physical con- 
ditions. Thus the lines arising in the cold component are 
narrow (Av— 1 km s _1 ) and the emission peak is located 
at the star position in FIRS 2. Besides, the kinetic tem- 
perature of this gas estimated from the NH3 lines is —13 
K in FIRS 2 and -28 K in LkHa 234. 

A warm envelope component is detected towards these 
sources traced by the emission of species like CS, CH3OH 
and H 2 CO. These species present enhanced abundances 
in regions where the icy grain mantles are evaporated 
(van der Tak et al~l l2000|) . They are also abundant in 
molecular ouflows where they can be released to the gas 



phase by shock fronts. In FIRS 2, the emission of the 
low energy transitions of these species arise mainly in the 
bipolar outflow. However, this association is not clear in 
LkHa 234 where their emisison could arise in the inner 
and warmer part of the envelope. 

Finally, we have strong evidences for the existence of 
a hot core in the Class protostar FIRS 2. The high den- 
sity measured at ambient velocities from the C S lines, 
the high temperature component of the CH3OH lines and, 
above all, the detection of the CH3CN lines with rotation 
temperature of —63 K show the existence of a hot core in 
this target. 

7.2. Physical and Chemical evolution of the YSOs 
envelopes 

7.2.1. Cold envelope 

In Table 9, we show the physical paramcnters and molec- 
ular column densities in the Class protostar FIRS 2 and 
LkHa 234. The molecules NH 3 , N 2 H+, H 13 CO+ and their 
deuterated compounds DCO + and N 2 D + trace the cold 
envelope component. The column densities of these species 
decrease by a factor of 5-10 from the Class protostar to 
the Type I Herbig Be star showing that the mass of the 
cold envelope decrease by at least a factor of 5 during the 
protostellar phase. Based on the NH3 and N 2 H + data we 
have also derived the kinetic temperature and size of this 
cold component. The kinetic temperature increases from 
—13 to —28 K and the size of the emitting region decreases 
from 21" to 6" from FIRS 2 a nd LkHa 234. This is consis- 
tent with previous results bv lFuente et al. ] ll2002l) which 
shows the protostellar envelope is dispersed and becomes 
warmer during its evolution to become a visible star. 

The warming of the cold envelope produces changes 
in its chemical composition. In Table 5 we show the 
N 2 H+/H 13 CO+ and NH 3 /N 2 H+ abundance ratios in sev- 
eral positions. Note that the N 2 H+/H 13 CO+ ratio is dif- 
ferent in the studied objects. The N 2 H + /H 13 CO + ratio is 
—17 in the cold young object FIRS 2 and —6 in the more 
evolved and warmer Herbig Be star LkHa 234. As dis- 
cussed in Section 3.2, this gradient in the N 2 H + /H 13 CO + 
ratio is very likely due to the H 13 CO + depletion in the cold 
envelope of the protostar. Molecular depletion is expected 
to be significant only for kinetic temperatures T^<20 K. 
Thus, molecular depletion is negligible in the envelope of 
LkHa 234 where the gas kinetic temperature is —28 K. 

Within the nitrogen chemistry, we have also studied 
the NH3/N 2 H + ratio. Since the beam is very different for 
the NH3 and N 2 H + observations, this abundance ratio is 
very dependent on the assumed source size. For LkHa 
234 we have been able to calculate both, the NH3 and 
N 2 H + emitting regions sizes. We have derived the same 
size, — 6"-8", for both molecules and the NH3/N 2 H + ra- 
tio is —10. In the case of FIRS 2 we have not been able 
to derive the source size from the NH3 emission but we 
have estimated a size of —21" from the N 2 H + observa- 
tions. We have assumed two limiting cases for the calcu- 
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Fig. 17. Histogram with the total beam-averaged molecular column densities estimated in the IM Class YSO NGC 
7129 - FIRS 2 and the HBe star LkHa 234. The column densities of the molecules tracing the cold envelope decrease 
at least by a factor of 5 between NGC 7129-FIRS 2 and LkHa234. An exception is H 13 CO+ that is depleted in the 
Class protostar (see text). The column densities of the molecules tracing the warm envelope and the PDR remain 
quite constant, in some cases increase, from the Class protostar to the HBe star. The SiO column density is lower 
by almost 2 orders of magnitude in LkHa 234 than in NGC 7129-FIRS 2. This is very likely related to the absence of 
an energetic bipolar outflow in the HBe star. 



lations of the NH3/N2H 4 " ratio in this source. Assuming 
a beam filling factor of ~1 for both molecules we obtain 
NHa/N2H + ^13 in this object. However if we assume that 
the size of the NH3 emission is ~ 21", like in the case of 
N2H+, we obtain NH3/N2H + ~33. In this case, we would 
have an NH3 abundance enhancement in the colder enve- 
lope of the Class protostar. Recent results in pre-stellar 
core show that the NH3 abundance could be enhanced in 
dense region s of these cores whe re the CO is expected to 
be depleted l|Tafalla et al. Il2004|) . 

7.2.2. Warm envelope 

When the star heats the envelope a sublimation front pro- 
ceeds outwards the star removing molecules from grain 
mantles. The region of the envelope in which the gas ki- 
netic temperature is high enough to evaporate the grain 
mantles is what we have called "warm envelope". The 



species released to gas phase are called "parent molecules" 
and increase significantly their abundances. The molecules 
CH3OH, NH3, and H2CO are within this group. These 
molecules drive a high temperature chemistry giving rise 
to "daughter" molecules like CH3CN. Within this scheme, 
CH3OH, H 2 CO and CH3CN are tracers of the warm 
part of the envelope where the ices have been evaporated 
l|Rodgers fc Charnlevl2 003) . Some of these species are also 
abundant in the molecular outflow where shock fronts re- 
move them from the grains mantles, and to a lesser ex- 
tent in the cold envelope. This is the case for NH3 with 
the emission of the low-lying transitions dominated by the 
cold envelope in both, FIRS 2 and LkHa 234, while the 
(3,3) line arise in the warmer component. 

In Table 9 we show the physical conditions and the 
molecular column densities in the warm envelope of 
FIRS 2 and LkHa 234. Contrary to the species tracing 
the cold envelope, the species tracing the warm enve- 
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lope present similar column densities in both targets (see 
Fig. 17). In fact, some species like H 2 CO and C 34 S seem 
to have larger column densities in LkHa 234. On the other 
hand, the size of the warm envelope is larger in the case 
of the Herbig Be star than in the case of the Class YSO. 
Thus, the mass and size of the warm envelope remains 
quite constant, or even increase, during the protostellar 
evolution. 

Although the column densities of the "warm envelope" 
species are not very different in the two YSOs, the ori- 
gin of their emission could be different. The CH 3 OH and 
H 2 CO emission seem to be dominated by the molecular 
in FIRS 2. In the case of CH 3 OH we have detected two 
components in FIRS 2. The one centered at ~7 km s _1 is 
clearly associated with the molecular outflow FIRS 2-out 1 
and dominates the CH3OH emission in all the transitions 
with E u <100 K. We have observed only two low-energy 
H2CO lines in our targets. The emission of these lines in 
FIRS 2 is clearly associated with the outflow FIRS 2-out 
1 (see Fig. 13). 

Intense CH3OH and H2CO lines have also been de- 
tected in LkHa 234. In this case the link between these 
species and the molecular outflow is not clear. In fact, 
these lines could arise from a warm inner envelope where 
the icy mantles are being evaporated releasing these 
species to the gas phase. Thus, although the CH3OH and 
H2CO column densities are very similar in both objects, 
the mechanism which removes these species from the icy 
mantles can be different in the Class protostar FIRS 2 
and the Herbig Be star LkHa 234. While shock fronts are 
very likely the main mechanism for the erosion of grain 
mantles in FIRS 2, ice evaporation could be a dominant 
mechanism in LkHa 234. 

7.2.3. CN,HCN 

The CN and HCN column densities increases from FIRS 2 
to LkHa 234 by a factor of ~ 6. This is easily understood 
within an evolutionary trend. Since the cold envelope has 
already been dispersed in LkHa 234, the UV radiation can 
penetrate deeper into the cloud affecting the chemistry all 
across the envelope in LkHa 234. The enhancement of 
ionization fraction (in particular, the enhancement in the 
C + abundance) in the LkHa 234 envelope could produce 
an enhancement in the fractional abundances of the nitro- 
genated chains CN and HCN relative to N 2 H+ and NH 3 . 

The CN/HCN ratio has been widely used as a 
PDR tracer thus far. Values of the CN/HCN ra- 
tio >1 has been considered as a prove of the exis- 
tence of a photon-domina t ed ch e mistry region (s e e e.g. 
Fuente et all Il993l fl995l l2003t iBachiller et all Il997t 
Rodriguez-Franco et al. 1 1998|) . We have detected a small 
region in the vicinity of the Class source in FIRS 2 in 
which the CN/HCN ratio is >1. This strongly suggests 
the existence of an incipient PDR in the inner part of this 
protostellar envelope. However, our limited angular reso- 
lution prevent us to discern if this PDR is associated with 



the Class source or with the IR companion star which 
is driving the outflow FIRS 2-out 2. 



7.2.4. Deuterated compounds 

We have observed DCO + and N2D" 1 " in order to be 
able to derive the D/H ratio in the cold envelope of 
these stars and its possible changes during the cen- 
tral object evolution. Enhancements of the D/H ratio 
over the 10 -5 ratio of [HD]/[H 2 ] have been found in 
dark clouds and young pr otostars (e.g. iButner et al. I 
1199.4 IWilliams et al~l Fl998V There are two main ways 
of producing these enhancements. Firstly, grai n surface 
chemi stry may enhanc e molecul ar D/H ratios |Tielensl 
EiH iRrown fcMillaJllflSflaN ICharnlev et al. I Il997h. 
Secondly, some key gas phase reactions involving destruc- 
tion of deuterated species run slower at low tempera- 
tures than the equivalent reactions with hydrogen, and 
this leads to molecular D/H enhancem ents where a cold 
gas phase chemistry has been active l|Roberts fe M illar 
l2000al iRoberts fe Millar"! l2000b|) . Furthermore, in colder 
gas, depiction of heavy molecules such as CO results 
m an increase of [H 2 D+1/[H+] and molecular D/H 
ratios foalgarno fc Leppl Il984t Brown fc Milkul Il989al 



ratios nu aigarno & t^ epp ii^»4t mrown &c Miiiarl iiM »aa: 
IRoberts fc Milla,^ l2000at IRoberts fc Millar I l2000hh. We 



have found DCO+/H 13 CO+-0.7 in FIRS 2 and -0.25 
in LkHa 234. Assuming a 12 C/ 13 C isotopic ratio of 89, 
this would imply DCO+/HCO+-0.008 in FIRS 2 and 
—0.003 in LkHa 234. These ratios are similar to that 
fo und in low mass Cla s s p rotostar IRAS+16293-2422 
bv lvan Dishoeck et al. I l|l995() b ut are lower that the val- 
ues found in dark clouds (see e.g. lTine et al. l200(il) . A cold 
gas chemistry as well as depletion of heavy molecules can 
explain the enhanced values DCO + /HCO + ratio found 
in the young stellar object FIRS 2. 

We have also observed the doubly deuterated 
formaldehyde D 2 CO. The lincwidths of the observed 
D 2 CO lines show that this species, like the non-deuterated 
compound H 2 CO, arise in the warm envelope. The D 2 CO 
column densities as well as the D 2 CO/H 2 CO ratio, are 
similar in LkHa 234 than in FIRS 2. This suggests that 
the deuterium fractionation remains constant (or even in- 
crease, see Table 9) in the warm envelope during the pro- 
tostellar evolution. Thus, the evolution of the deuterium 
fractionation is different in the cold and warm part of the 
envelope. Grain-surface chemistry may enhanced the deu- 
terium fractionation in the warm envelope where molec- 
ular evaporation is very likely the main chemical phc- 
nomcnum, while the cold gas chemistry and depletion 
could determine the evolution of the deuterium fractiona- 
tion in the cold envelope. 

The deuterium fractionation has been proposed as a 
chemical clock in YSOs. Our data confirms that the deu- 
terium fractionation changes significantly during the pro- 
tostellar evolution and consequently, can be used as a 
chemical clock. However, these changes are different in 
the cold and in the warm part of the envelope, and the 
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Table 9. Beam averaged column densities 





NGC 7129 - FIRS 2 


LkHa 234 


Age 


~ 3000 yr 


~ 10 5 yr 


n(cm~ 3 ) 


5 10 5 - 3 10 6 


2 10 5 




Cold envelope 




T fe 


13 K 


28 K 


n s 


~ 21 " 


~ 8" 


NH 3 


4.9 10 14 


4.0 10 13 


N 2 H+ 


3.8 10 13 


1.0 10 13 


H 13 CO+ 


2.2 10 12 


1.7 10 12 


DCO+ 


1.8 10 12 


4.2 10 11 


N 2 D+ 


5.4 10 11 


< 2.3 10 11 


N 2 H + /H 13 CO+ 


17 


6 


NH 3 /N 2 H+ 


13 -33 b 


10 


DCO + /H 13 CO + 


0.7 


0.2 


N 2 D+/N 2 H+ 


0.014 


< 0.02 




Warm envelope 




T fe 


> 50 K 


> 100 K 


n s 


~ 9" 


~ 14" 


c 34 s 


1.5 10 12 


1.7 10 12 


CH3OH 


9.8 10 14 


5.4 10 14 


ch 3 cn 


3.6 10 12 


> 2.0 10 12 


H 2 CO 


7.0 10 13 


8.0 10 13 


H 2 13 CO 


2.3 10 12 


2.4 10 12 


D 2 CO 


2.0 10 12 


3.7 10 12 


D 2 CO/H 2 CO 


0.03 


0.05 




PDR 




CN 


4.8 10 13 


6.1 10 13 


HCN 


1.6 10 13 


2.1 10 13 




Outflow 




SiO 


2.8 10 12 


2.0 10 11 


b Assuming beam filling factor ~ 1 and a source size of 2! 


the NH3 emitting 


region. 




Table 10. Chemical diagnostics for YSO evolution 




NGC 7129 - FIRS 2 


LkHa 234 


SiO/C 34 S 


- 2 


~ 0.1 


CN/N 2 H+ 


~ 1 


~ 6 


HCN/N 2 H+ 


- 0.4 


~ 2 


DCO+/HCO+ 


~ 0.008 


~ 0.002 


D 2 CO/DCO+ 


~ 1 


~10 



measured D/H ratio is very dependent on the molecular 
compounds used to determine it. 

8. Summary and Conclusions: Chemical-clocks in 
intermediate-mass YSOs 

We have carried out a molecular survey towards the IM 
YSOs FIRS 2 and LkHa 234. Our survey confirms that 
protostellar envelopes are very complex objects composed 
by several components characterized by different physical 
and chemical properties: 

1. A cold envelope observationally characterized by nar- 
row lines (Av~l km s^ 1 ) and low kinetic temperatures 



(Tfe < 50 K). This component is well traced by the 
species NH 3 ,N 2 H+, H 13 CO+, DCO+ and N 2 D+.' 

2. A warm envelope observationally characterized by 
large linewidths (Av>3 km s _1 ) and kinetic temper- 
atures (Tfc>50 K). The low energy transitions of CS, 
C 34 S, CH3OH, and H 2 CO arise in this warm compo- 
nent. 

3. A hot core characterized by high densities (n>10 6 
cm~ 3 ) and a high kinetic temperature (Tfe>100 k). 
The symmetric rotor CH3CN seems to be the best 
tracer of this hot component. 

4. In addition to these envelope components, energetic 
outflows are associated with very YSOs. In addition to 
CO and its isotopes, the SiO emission is very likely the 
best tracer of the outflow component for very young 
stellar objects. 

Once, we have used chemistry to determine the physi- 
cal structure of the YSOs, FIRS 2 and LkHa 234, we can 
determine the evolution of the protostellar envelopes of 
IM stars during the protostellar phase. FIRS 2 is an IM 
Class objects while LkHa 234 is a very young (and still 
deeply embedded) HBe star. As expected, different phys- 
ical conditions and chemistry are found in these objects. 
The Class IM is a cold object (T fc - 13 K) in which 
molecular depletion is still important. We have no evi- 
dence of a change in the N 2 H+ abundance between both 
objects. The decrease in the N 2 H + column density ob- 
served in Fig. 17 is the consequence of the disruption of 
the cold envelope during the stellar evolution. However, 
we have found evidences for H 13 CO + depletion in FIRS 2. 
Molecules like CH 3 OH and H 2 CO are expected to trace 
mainly the warm region in which grain mantles have been 
evaporated. The column densities of these molecules re- 
main constant (or increase) from the Class IM to the 
more evolved HBe star, suggesting that the abundances 
of these molecules and/or the mass of the warm gas in- 
creases with the protostellar evolution. The detection of 
CH 3 CN and the high temperature derived from it (Tfc~63 
K) shows that a hot core has developed in the Class pro- 
tostar FIRS 2. Thus, our results suggest an evolutionary 
sequence in which as the protostar evolves to become a 
visible star, the total column density of gas is decreas- 
ing while the amount of warm gas remains quite constant 
or slightly increase. These physical and chemical changes 
imply important changes in the beam averaged column 
densities during the protostellar evolution. 

Based on our observational study of FIRS 2 and LkHa 
234, we propose some abundance ratios that can be used 
as chemical clocks in YSOs. These ratios are defined to 
be useful tools to discern between different evolutionary 
stages of YSOs, but do not correspond to the actual abun- 
dance ratios in any of the envelope components. They have 
been calculated with beam averaged column densities and, 
as largely discussed in this paper, are the consequence of 
a complex physical and chemical evolution in the whole 
envelope. 
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In Tabic 10 we list the proposed chemical clocks. We 
find the maximum variation between the Class and the 
HBe object when comparing the SiO/C 34 S ratio. This ra- 
tio decreases by a factor of ~ 20 between these two objects. 
This decrease is the consequence of the decay of the bipo- 
lar outflow phenomenum during the protostellar evolution. 
This ratio is specially useful to determine the evolutionary 
stage of the youngest objects which are associated with 
the most energetic bipolar outflows. The nitrogen chem- 
istry is also useful to determine the evolutionary stage of 
YSOs. The CN/N 2 H+ and HCN/N 2 H+ ratios are larger 
by a factor of <~ 6 in the HBe star than in the Class ob- 
ject. As commented above this is mainly due to the fact 
that the LkHa234 envelope is thinner and warmer than 
that of the FIRS 2. These ratios would probably be more 
useful to differentiate between objects in the late proto- 
stellar evolution when the protostellar envelope becomes 
optically thinner. Finally, the deuterated species could 
also be a good indicator of the protostellar evolution. The 
DCO + /HCO + ratio decreases by a factor ~ 4 due to the 
warmer envelope in LkHa 234. However, we can have a dif- 
ferent behavior in the deuterated species whose emission 
arise in the warm envelope component. This is the case of 
the doubly deuterated compound D2CO. We obtain the 
largest contrast in the abundance ratio if we compare the 
D 2 CO and DCO+ abundances. The D 2 CO/DCO+ ratio 
increases by a factor of 10 from the Class to the Herbig 
Be star. However, we should be cautious using this ra- 
tio because we are comparing species arising in different 
regions of the envelope. 
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Fig. 8. In the central panel we show the integrated line 
intensity map of the SiO J=2^1 line toward FIRS 2. 
Contour levels are from 0.25 K kms -1 to 4 K kms -1 by 
steps of 0.5 K kms -1 . Surrounding the central panel, we 
show the observed spectra of the 13 CO 1— »0, SiO 2— >1 and 
H 13 CO + 1— >0 lines towards selected positions of outflows 
FIRS 2-out 1 (top panels) and FIRS 2-out 2 (bottom pan- 
els). The intensity scale in the spectra is main brightness 
temperature. 



A. Fuente et al.: Chemical evolution in the environment of intermediate mass young stellar objects: 



23 




Fig. 10. Observed spectra of CH3OH and H 2 CO lines 
towards the (0,0) position in NGC 7129 FIRS 2 and 
LkHa 234.. 
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Fig. 11. The (a) panel shows the rotational diagram for 
CH 3 OH in the (0,0) position of NGC 7129 - FIRS 2. 
We have degraded the angular resolution of the 1.3mm 
maps (the 5^4 and 8^7 methanol lines) in order to have 
the angular resolution of the 3mm maps (lines 2— >1 and 
11 — >10 lines). Two rotation temperatures (T rot =17 K and 
T rot >80 K are required to fit the data). In the (b) and (c) 
panels we plot the line velocitcs and lincwidths, V and AV, 
of the CH3OH lines as given by the gaussian fits shown 
in Table 3 versus the upper state energy of the observed 
transition (E M ). Finally, the (d) panel shows V versus Av 
for the CH3OH lines. This panel clearly shows the exis- 
tence of two veloctiy components, one at ^—10 kms^ 1 and 
the other at ^—6 kms^ 1 in this position. 



Fig. 12. The top panels show the rotational diagram for 
CH3OH in LkHa 234 in the two limit cases assumed in 
the paper, panel (a) corresponds to a beam filling factor 
~ land panel (b) corresponds to a point source. In the 
bottom panels we plot the linewidths, Av, derived from 
the Gaussian fits versus the upper energy state of the tran- 
sition, E u , and the central velocity, V, for all the lines. In 
the case of LkHa 234 we have only one velocity compo- 
nent. 
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Fig. 14. Rotational diagram of CH3CN in the (0,0) posi- 
tion of FIRS 2 and LkHa 234. 
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Fig. 15. Observed spectra towards the (0,0) position in 
LkHa 234 and NGC 7129 - FIRS 2. The warming of the 
gas during the protostellar evolution, produces significant 
chemical changes specially in the cold component. 
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Fig. 16. CN 1— >0 /HCN l->0 integrated intensity ratio in 
NGC 7129 - FIRS 2 (left) and LkHa 234 (right). Contour 
levels are 0.10 to 0.35 by 0.05 in NGC 7129 - FIRS 2 and 
0.35 to 1.1 by 0.05 in LkHa 234. 



